Disruption of reflecting Bose-Einstein condensates due to inter-atomic interactions 

and quantum noise 
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We perform fully three-dimensional simulations, using the truncated Wigner method, to inves- 
tigate the reflection of Bose-Einstein condensates from abrupt potential barriers. We show that 
the inter-atomic interactions can disrupt the internal structure of a cigar-shaped cloud with a high 
atom density at low approach velocities, damping the center-of-mass motion and generating vor- 
tices. Furthermore, by incorporating quantum noise we show that scattering halos form at high 
approach velocities, causing an associated condensate depletion. We compare our results to recent 
experimental observations. 
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I. INTRODUCTION 

Recent experiments P, |3l have demonstrated quantum 
reflection of Bose-Einstein condensates (BECs) from a 
silicon surface, and showed the potential applications of 
semiconductor surfaces as atom mirrors and traps. This 
is very exciting in the light of the intense interest inpre- 
cise manipulation of cold atoms with atom chips |3, ll| , 
which have been used to reflect, transport, split, and re- 
combine atom clouds H H II |8, J., ifl]. In principle, 
semiconductor surfaces have many advantages over con- 
ventional atom chips: they can be precisely patterned, 
require no magnetic fields or currents, and hence avoid 
the problems of fragmentation • In order to optomize 
the effective manipulation of atoms, theoretical work is 
required to understand the reflection process, and in par- 
ticular, how it may disrupt the cloud. Previous theoreti- 
cal work [T^ , based on solving the cylindrically symmet- 
ric Gross-Pitaevskii equation (GPE), showed that inter- 
atomic interactions can disrupt the internal structure 
of the cloud on reflection, creating topological excita- 
tions such as solitons and vortices. Recent theoretical 
work on colliding BECs Ts', however, demonstrated 
that the bare GPE is insufficient for describing the pro- 
duction of scattering halos and the associated depletion 
of the BEG, develo ping the truncated Wigner method 
(TWM) [H El El llal , which includes quantum noise, 
as an alternative method to model this process. Further- 
more, the experimental results have confirmed that large 
scattering halos are indeed created for particular param- 
eter regimes 0- 

In this paper, we reproduce the experimental results 
by performing fully three-dimensional simulations using 
the TWM. In contrast to the bare GPE, this method 
models quantum fluctuations of the condensate, and can 
hence describe processes which require spontaneous ini- 
tiation. One such process is the pairwise scattering of 
condensate atoms into unoccupied modes, which is the 
precursor for the formation of a scattering halo. Once 
these previously unoccupied modes have some finite oc- 



cupation, further atoms enter them via Bosonic stimula- 
tion, creating a scattering halo. By comparing multiple 
simulations the TWM may also be used to calculate the 
condensate depletion associated with the halo. 

We explore the role of the cloud geometry, density and 
approach speed in the reflection of BEGs from abrupt 
potential steps. We find that the disruption due to the 
inter- atomic interactions, which we refer to as interfer- 
ential disruption, is most pronounced for cigar-shaped 
BEGs with high atom densities, approaching the poten- 
tial barrier along its long axis at low velocities. This ef- 
fect generates vortices and causes an associated damping 
of the center-of-mass motion. The production of scatter- 
ing halos is also enhanced by elongating the cloud and 
increasing the density, but in contrast to the interfer- 
ential disruption, is most pronounced at high approach 
velocity. As mentioned previously, the creation of the 
scattering halo causes a depletion of the condensate. We 
compare our theoretical results to recent experimental 
observations of excited reflected clouds and scattering ha- 
los Hi. 

We study two sets of ^'^Na BEG parameters, referred to 
as A and B, taken from recent experiments on quantum 
reflection of BEGs from silicon surfaces. BEG B is more 
than twice as dense (its equilibrium peak density no is 
5.2 X 10^2 cm-3) as BEG A (no = 2.0 x lO^^ cm'^). Both 
BEGs have different widths in each coordinate direction. 
The long axis of BEG B is always perpendicular to the 
barrier, but we examine the dynamics of BEG A in three 
orientations, in which the cloud is accelerated in each of 
the three principal coordinate directions. By studying 
the behaviour of a single BEG in different orientations 
we draw conclusions about the role of BEG shape and 
orientation in determining the subsequent dynamics. We 
show that our theoretical predictions for the two sets of 
parameters are in good agreement with the experimental 
observations. 

The paper is organised as follows: in section ^ we 
discuss the theoretical model used to simulate the BEG 
dynamics, in section UlTl we present results for BEG A, in 
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section llVl we present results for BEC B, and in sectionlVI 
we summarize our findings and conclude. 



II. THEORETICAL MODEL OF THE BEC 
DYNAMICS AND QUANTUM FLUCTUATIONS 
OF THE FIELD 

A. The truncated Wigner method 

Heuristically, the TWM can be thought of as a clas- 
sical field technique which simulates quantum vacuum 
fluctuations by adding appropriate classical random fluc- 
tuations to the coherent field of the EEC's intial state. 
The added fluctuations are referred to as virtual particles, 
and the atoms in the BEC initial state are referred to as 
the real particles. This approach is valid for high densities 
of real particles. Specifically, in a homogenous system, 
the number of real particles must be much greater than 
the number of virtual particles. The derivation of this 
method is outlined in Ref. |0| and described in detail 
in Ref. 0|. The dynamical equations of the TWM are 
identical to the projected GPE \it\ . 

The wavefunction (x, t) is modelled by the mode ex- 
pansion 



j=M 



* (x, = ^ aj {t) e 



(1) 



where t is time, V — L^LyLz is the volume contained 
in the coordinate space and aj (t) is the amplitude of 
the mode with wavevector kj, normalised such that 



J2jZi^ a* (t) OLj [t) is the total number of atoms IS! . The 
mode space is spherical, being contained within a maxi- 
mum cut-off wavevector modulus to prevent Fourier alias- 
ing 1^^. The wavevectors are defined as 



in which a = 2.9 nm is the s-wave scattering length, using 
an imaginary time algorithm |19| | . The quantum fluctua- 
tions are introduced by combining this real particle field 
ip (x) with a field of virtual particles x (x) to create the 
total field ^ (x, 0) . The virtual particle field is defined as 



X(x) 



7f E 



(6) 



in which the complex amplitudes Xj have a Gaussian dis- 
tribution with the properties iXiXj) = 5% ^-iid iXiXj) — 
0, Q| . This means that on average each mode is pop- 
ulated by half a virtual particle, so that the total number 
of virtual particles is approximately M/2. 

The derivation in Ref. jT3| yields the stochastic dif- 
ferential equation for each mode within the low energy 
subspace 



1 ■' 



1 



dt 2m yjy 

(7) 

where VcyA is the total external potential. We compute 
the dynamics of the BEC by solving these equations using 
the Fourth- order Runqe-Kutta in the Interaction Picture 
algorithm (RK41P) 20]. 



B. Application to the reflection problem 

At t — the magnetic trap is displaced by a distance 
Aa; along the x-direction, such that it is now centered 
at (0, 0, 0), creating the new potential profile ?7trap- This 
accelerates the BEC towards an abrupt potential step 
C^waii of height V , which is positive, in the y — z plane at 
X = 0, so that 



Lx Ly 



(2) 



where Pj, Qj and rj are integers. 

For t < 0, the BEC is held in a magnetic trap of fre- 
quencies ujx, ujy and ujz, centered at (—Ax, 0,0), with 
potential profile 



'-'tiap — 2 



(3) 



where m is the mass of a single ^'^Na atom. The initial 
state is determined by solving the three-dimensional GPE 



fc2y2 

-^— + Utr.p + Uo\^\^ 

2m 



where 



Uo 



m 



(4) 



(5) 



V x>0 
a; < 



(8) 



Hence, the impact velocity of the cloud at the wall Vx ~ 
LJxAx j23|. If the mean kinetic energy of the atoms at 
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the barrier (E) « <C V, all atoms are reflected. If 
(E) > V, there is finite transmission. 

In this paper we compare our theoretical predictions 
to experimental obvservations of quantum reflection from 
the Casimir- Polder potential 0, |2| . Unfortunately, it is 
impossible to model the Casimir-Polder potential directly 
as this accelerates atoms to large wavevectors which ex- 
ceed our maximum cut-off wavevector, which is limited 
by computer memory. However, since the Casimir-Polder 
potential, and particularly that of a pillared surface, 
varies over a distance of less than the healing length , 
it can be approximated to an abrupt potential barrier. 
Furthermore, previous theoretical work [T^ has demon- 
strated that the BEC dynamics are qualitatively inde- 
pendent of the form of the reflecting potential barrier. 



1. Absorption of transmitted real particles 

Transmitted atoms are absorbed by an imaginary po- 
tential, given by 



^imag — 



-Cx X > Q 
a; < 



(9) 



where C is a positive constant. The imaginary potential 
introduces a damping term into the equations for the 
mode amplitudes, which removes atoms for x > Q. How- 
ever, the imaginary potential must only absorb the real 
particles, not the virtual particles, hence maintaining the 
appropriate amount of quantum noise. We overcome this 
problem by calculating the damping rate of the imaginary 
potential, and adding extra virtual particles by increas- 
ing each mode amplitude Uj by 



(10) 



every time step At j22|, where /S.V is the volume around 
one spatial grid point, and the complex spatial process 
has the properties {W^) = 0, (W^W^) = 0, and 

Hence the total external potential is t/ext = Ui 

^wall H~ f^imag- 



III. DYNAMICS OF BEC A 
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For BEC A, N ^3x10^ 



— 2Tr X 3.3 rad s 



27T X 6.5 rad s 
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For these 



27r X 2.5 rad s~ , and w_ 

parameters uq — 2.0 x IQ-'^^ cm~'^. FigureCJa) shows two 
views of a constant density surface of the BEC initial 
state at i = without added quantum fluctuations (the 
axes are included in the figure). In this section V = 10~^^ 

2 

J, which far exceeds ^-^^ for all considered, and hence 
all atoms are reflected. 

Since the cloud is not symmetric, we can explore the 
role of the BEC geometry by accelerating it in different 
directions. The parameters as defined above will be re- 
ferred to as orientation 1. We will later rotate the cloud 



such that ujx — 2t: X 2.5 rad s 



2tt X 3.3 rad s 



and ujz = 2Tr X 6.5 rad s"^, hence presenting its long axis 
towards the potential barrier. This will be referred to 
as orientation 2. Finally, we will rotate the cloud such 



that 



27r X 6.5 rad s 



27r X 2.5 rad s , and 



= 27r X 3.3 rad s"""^, hence presenting its short axis 
towards the potential barrier. This will be referred to as 
orientation 3. 



A. Interferential disruption 

We initially perform simulations without added quan- 
tum noise in order to study the interferential disruption. 




FIG. 1: Constant density surfaces of BEC A without quantum 
fluctuations in orientation 1, reflecting from a barrier of height 
V = 10"^° J for = 1.2 mm s"^ at i = (a), 75 (b), 100 (c) 
and 125 ms (d). The left and right-hand columns shows views 
from and towards the barrier respectively [axes are shown in 
(a)]. The vertical bar shows the scale. 



Figure n shows constant density surfaces for BEC A in 
orientation 1 without quantum fluctuations undergoing 
reflection from the potential barrier at i = 75 (b), 100 (c), 
and 125 ms (d), for — 1.2 mms^^. The left and right- 
hand columns show views of the cloud from and towards 
the wall respectively [see axes in Fig.^a)]. Figure ^b) 
shows the mid-point of the oscillation, at which time the 
cloud is modulated by the standing wave created by the 
superposition of the incident and reflected matter waves. 
The left-hand flgure shows the flat constant density sur- 
face against the potential barrier. Previous theoretical 
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work |l2l | demonstrated that the inter-atomic interactions 
in the high-density peaks of the standing wave cause the 
momentum of some atoms to be transferred into the ra- 
dial direction. These atoms appeared as jet-hke "lobes" . 
In three dimensions, the "lobes" can be seen to form more 
of a "doughnut" [Fig.njc)]. The trap causes the "dough- 
nut" to collapse as it rebounds away from the wall, dis- 
rupting the internal structure of the cloud [Fig. Did)]. 
At the end of the oscillation {t = 150 ms w tt/lOx), the 
radial momentum of the "doughnut" has created vortex 
lines which pierce the BEC, roughly parallel to the j/-axis 
[Fig.|2Ia)]. The circular arrow in Fig.|21[a) shows the di- 
rection of quantized circulation. These theoretical find- 
ings are in agreement with experimental observations of 
"excited and sometimes fragmented" reflected clouds Q . 

The excitation of the internal structure of the cloud 
removes energy from the longitudinal motion, and hence 
this effect is associated with a damping of the center-of- 
mass motion. The upper solid line in Fig. EJg) is the 
mean BEC x position (x) as a function of t for these 
parameters, and it is clearly strongly damped. The effect 
is described as interferential disruption, since it arises 
from an interference pattern, in this case produced by 
the superposition of the incident and reflected matter 
waves. 

The "lobes" form if atoms can move radially outwards 
from the center of the BEC to its extremities before the 
reflection process is complete The time taken for 

this to occur has the lower limit ti = Ir/vs, where Ir is 
the radial half-width of the cloud and Vs ~ yvi^noo/Trrr? 
is the average speed of sound, if the atoms are to travel 
through the cloud without exciting the BEC. The reflec- 
tion time is t^ « Ix/vx, where Ix is the BEC's longitudi- 
nal half-width. Consequently we expect to observe lobes 
and subsequent disruption if 

Vx <Vs{lx/lr) ■ (11) 

For the simulation described above, we have the in- 
teresting situation that equation 1111 is satisfied if Ir is 
taken to be Iz, but not if Ir is taken to be ly. Con- 
sequently, the "doughnut" is particularly pronounced in 
the z-direction, and hence its momentum is largest in the 
^-direction. This explains why the resulting vortex lines 
are roughly parallel to the y-axis. 

As Vx increases equation ^2 is no longer satisfied. Fig- 
ure [2Ib) shows the reflected cloud in orientation 1 at 
t — 150 ms for Vx — 2.1 mm s^^. For these parame- 
ters the "doughnut" does not form. Hence, at the end 
of the oscillation the cloud contains no topological exci- 
tations and its appearance is similar to the initial state 
at t = 0. Consequently, there is little damping of the 
center-of-mass motion, as shown by the lower solid curve 
in Fig.Hg). 

We may control the disruption by changing the geom- 
etry of the cloud, i.e. by changing the ratio Ix/lr in equa- 
tion ^2 The ratio Ix/lr is largest for orientation 2, and 
smallest for orientation 3, so we would expect more and 




t (ms) 200 

FIG. 2: Constant density surfaces of BEC A without quantum 
fluctuations in orientation 1 [(a) and (b)], orientation 2 [(c) 
and (d)], and orientation 3 [(e) and (f)] after one complete os- 
cillation, having reflected from a barrier of height V — 10~^° 
J. The left and right-hand figures show results for Vx = 1.2 
and 2.1 mm s~^ respectively. The circular arrow in Fig. (a) in- 
dicates the direction of circulation around a vortex line. The 
vertical bar shows the scale, the axes are shown in (a), (g); 
(x) versus t for orientation 1 (solid curves), orientation 2 (dot- 
ted curves) and orientation 3 (dashed curves). In all cases the 
upper and lower curves correspond to Vx = 1.2 and 2.1 mm 
s^^ respectively. 



less disruption for these parameters respectively. This 
is supported by our numerical simulations. The BEC in 
orientation 2 is highly disrupted for Vx = 1.2 mm s~^ and 
contains many topological excitations [Fig. [2Ic)] , and is 
mildly disrupted for Vx — 2.1 mm s^^ [Fig.|2Id)]. How- 
ever, the BEC in orientation 3 shows virtually no disrup- 
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tion at either impact velocity [Fig.jSJe) and (f)]. Again 
we find that severe disruption of the BEC internal struc- 
ture is associated with strong damping of the center-of- 
mass motion: the damping is largest for orientation 2 
[dotted curves in Fig. [2Ig)], and smallest for orientation 
3 [dashed curves in Fig. Etg)] • These results illustrate 
that we expect severe interferential disruption for large 
Ix and small Ir, i.e. a cigar-shaped BEC with its long axis 
perpendicular to the barrier, as predicted by eauation llll 
Conversely, we expect mild interferential disruption for 
small Lj; and large Z^, i-e. a pancake-shaped BEC with its 
short axis perpendicular to the barrier. 



B. Inclusion of quantum fluctuations 

We also investigated the role of quantum fluctuations 
of the field by adding noise to the initial state at t = 0, 
as described in section ^ Figure 13 shows the resulting 
simulated absorption images of the BEC in the y — x 
plane after one complete oscillation for various orienta- 
tions and approach velocities. The interferential disrup- 
tion has least effect in this plane for the reasons explained 
in section IlII Al For the parameters taken from the first 
BEC quantum refiection experiments [ij (orientation 1), 
the scattering halo is very weak for an approach velocity 
of 1.2 mm [Fig. OJa)]. This is unsurprising as scat- 
tering halos were not observed in the experiments. For 
the larger approach velocity of 2.1 mm s^^ [Fig. Olb)] 
the halo is much clearer, but still relatively faint. This 
result is confirmed in Fig. ^g) , which shows the coherent 
number of atoms Nc as function of t for Vx = 1.2 and 
2.1 mm s~^ (upper and lower solid curves respectively), 
calculated by the method described in Ref. [lj|. As the 
BEC refiects, coherent atoms from the BEC are trans- 
ferred into the incoherent scattering halo, depleting the 
condensate. The final condensate fraction is much less 
(57% compared to 88%) for the larger of 2.1 mm s~^ 
than for — 1.2 mm s^^. 

If the BEC is orientated such that its long axis is per- 
pendicular to the wall (orientation 2), the scattering halo 
is much more pronounced [Fig.|3fc) and (d)], and there is 
a greater associated depletion of the condensate [dotted 
lines in Fig.|2Ig)]. Conversely, the halo is only just visible 
for the BEC in orientation 3 [Fig. |2Ie) and (f)]. These 
results indicate that the production of scattering halos 
is most pronounced for cigar-shaped BECs, approaching 
the potential barrier along its long axis with large Vx- 
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FIG. 3: Simulated absorption images of BEC A including 
quantum fluctuations in orientation 1 [(a) and (b)], orienta- 
tion 2 [(c) and (d)], and orientation 3 [(e) and (f)] after one 
complete oscillation, having reflected from a barrier of height 
V = 10"'^'' J. The left and right-hand flgures show results 
for Vx ~ 1.2 and 2.1 mm s~^ respectively. The horizontal 
bar shows the scale, and the dashed vertical lines indicate the 
position of the barrier, (g): Nc versus t for orientation 1 
(solid curves), orientation 2 (dotted curves) and orientation 
3 (dashed curves). In all cases the upper and lower curves 
correspond to Vx = 1-2 and 2.1 mm s~^ respectively. 



For BEC B, N ^ 10^ tOx ^ ^tt x 4.2 rad s~\ 
ujy — 2tt X 5.0 rad s^^, and ujz — 2tt x 8.2 rad s~^. These 
parameters give no = 5.2 x 10^^ cm~^, which is more 
than double that in BEC A. Figure ^a.) shows a con- 
stant density surface of BEC B without added quantum 
fiuctuations at t = 0. In this section, V is reduced below 



^^^^ , in order to investigate the effect of finite transmis- 
sion. Transmitted atoms are absorbed by an imaginary 
potential, as described in section ITll 
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FIG. 4: Constant density surfaces of BEC B without quantum 
fluctuations reflecting from a barrier of height V = 7.97 x 
10"^^ J for = 2.1 mm s"^ at f = (a), 70 (b) and 120 ms 
(d). The axes are shown in (a), (c): Shce through the BEC in 
the y — X plane at t = 70 ms. The arrows show the direction 
of circulation around a vortex ring, (e): Simulated absorption 
image of BEC B including quantum fluctuations in the y — x 
plane at t = 120 ms. The arrow in Fig. (a) indicates the scale. 



A. Low approach velocity 

Initially, we set V = 7.97 x lO'^^ J and = 2.1 
mins"^, for which we would expect a reflection probabiliy 
of 0.5 for a non-interacting plane wave. 

Since the atom density is higher in BEC B, and its long 
axis is perpendicular to the barrier, equation II II predicts 
more pronounced interferential disruption. This is con- 
firmed by the numerical simulations. Figures Elb) and 



(d) show constant density surfaces of BEC B without 
quantum fluctuations undergoing reflection from the po- 
tential barrier at i = 70 and 120 ms respectively. At 
t = 70 ms the standing wave has formed, but the inter- 
atomic interactions have already caused disruption of the 
internal structure and "doughnut" formation. In this 
high-density regime, some vortices are formed directly 
from the modulation of the standing wave, as has been 
demonstrated in optical lattices The slice through 
the y — X plane in Fig. QJc) shows the bending of the 
nodal lines, which is characteristic of the snake insta- 
bility 0ili3> before the "doughnut" has collapsed. A 
vortex ring has formed where one nodal line is bending 
and breaking, and is encircled by two arrows, indicat- 
ing the direction of quantized circulation. At the end of 
the oscillation [Fig. El^d)], the internal structure of the 
cloud is highly disrupted, and contains many topologi- 
cal exciations. This is in agreement with the results of 
Ref. 0, which reports the observation of disruption to 
the BEC internal structure for approach velocities of ap- 
proximately 2 mm s~^ and below. 

A corresponding simulation including quantum fluctu- 
ations shows that the higher density also causes more 
pronounced halo formation. Figure EJe) is a simulated 
absorbtion image of the cloud in the y—x plane a,tt— 120 
ms, showing a very clear, large and dense halo. 

We obtain qualitatively similar results in the case of 
zero transmission, by setting V — 10^"^° J, as in sec- 
tion |llll 



B. High approach velocity and comparison with 
experiment 

Finally, we increase V to 1.67 x 10~^^ and Vx to 3 mm 
s~^ in order to compare our predictions directly with 
experimental results [2] . For these parameters the plane- 
wave reflection probability is 0.5, approximately the same 
as that observed in experiment at this on reflection 
from the Casimir-Polder potential of a pillared silicon 
surface. 

At this larger velocity there is very little interferential 
disruption, so without quantum fluctuations the reflected 
cloud has a reasonably smooth density profile (Fig. |31 
upper inset). However, when we include quantum fluc- 
tuations an even larger scattering halo is produced, as 
shown in Fig.|31 lower inset. As before, we obtain qual- 
itatively similar results in the case of zero transmission 
for V = 10-30 J. 

Our simulations agree with the results of Ref. 2] , which 
reports the observation of large scattering halos for an 
approach velocity of 3 mm s"^, but no interferential dis- 
ruption for > 2 mms"^. An experimental image is 
shown in Fig. |S1 for comparison with our simulated ab- 
sorption image. The experimental and theoretical im- 
ages both show a dense and coherent cloud that has 
reflected cleanly from the barrier, superimposed on the 
background of the incoherent and comparitively dilute 
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FIG. 5: Experimental absorption image of BEG Bfor Vx = 3.0 
mm at t = 120 ms, having reflected from the Gasimir- 
Polder potential of a pillared silicon surface. The field of view 
is 500 fj,m, the vertical dashed line indicates the position of 
the barrier. Lower inset: corresponding simulated absorption 
image in the y — x plane including quantum fluctuations for 
refiection from a barrier of height V = 1.67 x 10^'^^ J. Upper 
inset: equivalent constant density surface excluding quantum 
fluctuations, axes are shown in the figure. 

halo. These findings demonstrate that the production of 
scattering hales is most pronounced for dense BECs with 
large v^. 



V. CONCLUSIONS 



It has already been demonstrated that abrupt poten- 
tial barriers, such as the Casimir-Polder potential of a sil- 
icon surface, could be exploited in applications as atom 
mirrors and traps [H. Furthermore, surface potentials 
can be routinely patterned in order to create devices to 
manipulate BECs in more imaginative ways. The results 
of this paper are more generally relevant to any reflection 
process, for example using magnetic mirrors lal, sheets of 
laser light, or arrays of current-carrying wires |lOj |. These 
techniques can also be extended to create more exotic po- 
tentials, such as diffraction gratings [2^l27l |. 

Our findings suggest that these experiments may be 
most successful for dilute, pancake-shaped BECs, with 
their short axis perpendicular to the barrier. Further- 
more, our results show that there are two competing ef- 
fects which may spoil the smooth profiles of the incident 
BECs and produce excitations: interferential disruption 
and the production of scattering halos. The former is 
most pronounced at low approach velocity, the latter at 
high approach velocity. Applications of BEC reflection 
may be most effective for moderate approach velocities 
which are above the threshold for interferential disrup- 
tion, but below the threshold for the formation of a scat- 
tering halo. Any future devices based on reflection must 
also account for possible damping of the center-of-mass 
oscillations associated with the production of excitations. 

We thank T. Pasquini for helpfifl discussions, and his 
kind permission to include the experimental image. This 
work was the supported by the Royal Society (London), 
and the Marsden Fund under contracts UOO-0590 and 
UOO-323. 
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